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in pp Collisions 
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Abstract 

We present a search for electroweak production of single top quarks in 
the electron+jets and muon+jets decay channels. The measurements use 
~ 90 pb~^ of data from Run 1 of the Fermilab Tevatron collider, collected 
at 1.8 TeV with the D0 detector between 1992 and 1995. We use events 
that include a tagging muon, implying the presence of a 6 jet, to set an upper 
limit at the 95% confidence level on the cross section for the s-channel process 
pp tb + X of 39 pb. The upper limit for the t-channel process pp tqb + X 
is 58 pb. 



*PACS numbers: 14.65.Ha, 12.15.Ji, 13.85. Qk 
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The top quark is the charge +2/3 weak-isospin partner of the bottom quark in 
the third generation of fermions of the standard model (SM). It is extremely massive 
at 174.3 ±5.1 GeV and, with an expected width of 1.5 GeV [0], it decays before 
hadronization almost exclusively into a W boson and a b quark. At the Tevatron pp collider, 
most top quarks are pair-produced via the strong interaction through an intermediate gluon. 
This was the mode used in its observation [Q and subsequent studies of its properties, 
including measurements of the ti production cross section of 5.9 ± 1.7 pb by the D0 collab- 
oration 0, and 6.5tJ;J pb by the CDF collaboration 0. A second production mode is 
predicted to exist, where top quarks are created singly through an electroweak Wtb vertex [Q. 
Many processes beyond the SM can boost the single top quark cross section [|^. In the 
absense of a cross section excess, measurement of the electroweak production of single top 
quarks could provide the magnitude of the CKM matrix element Vtb [§[, since the cross 
section is proportional to |V^tf,|^. In this Letter, we describe a search for single top quarks at 
the Tevatron, using data collected from 1992-1995 at app center-of-mass energy of 1.8 TeV. 



u d u d 




FIG. 1. Leading order Feynman diagrams for single top quark production at the 
Tevatron, where (a) shows the s-channel mode, and (b) the t-channel mode. 



The standard model predicts three modes for the production of single top quarks at a 
hadron collider. The first is the s-channel process q'q^tb, illustrated in Fig. |T](a). For a top 
quark mass rrit of 175 GeV, this has a cross section calculated at next-to-leading-order (NLO) 
of 0.73 ± 0.10 pb Following the decay of the top quark, these events contain a W boson 
and two b quarks that hadronize into two central jets with high transverse momentum (pt)- 
The second production mode, shown in Fig. ^(b) and sometimes referred to as 'W-gluon 



fusion", is a t-channel process, q'g—ytqb. The NLO cross section is 1.70 ± 0.24 pb IT^. This 
process produces a W boson, a forward light-quark jet, and two central b jets, one with 
high pt and the other with low pt- We have searched for both production modes, with 
decay of the W boson into eu or fiu, and identification of a 6 jet via a tagging muon. A third 
mode occurs via both the s-channel and t-channel, bg^tW, with a final state containing 
two W bosons and a single b jet. The leading-order cross section for this process is only 
0.15 pb [0, and, with ^ 90 pb~^ of available data, there is no possibility of separating it 
from background. Throughout this Letter, we use Hb" to refer to both tb and the charge- 
conjugate process ib, and to both tqb and iqb. 



The D0 detector [|T^ has three major components: a central tracking system including 
a transition radiation detector (TRD), a uranium/liquid- argon calorimeter, and a muon 
spectrometer. For the measurement in the electron channel, we use 91.9 ±4.1 pb~^ of 
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data collected with a trigger that required an electromagnetic (EM) energy cluster in the 
calorimeter, a jet, and missing transverse momentum (-^t)- For events passing the final 
selection, the efficiency of the trigger is 90-93%, depending on the location of the EM 
cluster in the calorimeter. In the muon channel, we use 88.0 ± 3.9 pb~^ of data acquired 
with several triggers, which required or a muon with a jet. The combined efficiency of 
these triggers is 96-99%. A third data sample, obtained with a trigger requiring just three 
jets, is used for measuring one of the backgrounds. Since the multijet cross section is very 
large, this trigger was prescaled, and we have 0.8 pb~^ of such data. Each of the three 
samples contains approximately one million events. 

To determine whether an EM energy cluster was generated by an electron, we require it to 
be isolated from other activity in the calorimeter and use a five- variable likelihood function 
to discriminate electrons from background. This likelihood includes the fraction of cluster 
energy contained in the EM region of the calorimeter (> 90% for electrons), the cluster shape 
(it must resemble an electron and not a pion), the presence of a well-matched track between 
the cluster and a primary interaction vertex (to discriminate against photons), the dE/dx 
energy loss along the track (consistent with a single particle and not from a photon conversion 
into a pair of charged particles), and the TRD response (matching that of an electron and 
not a pion). An electron is then required to have transverse energy Et > 20 GeV, and to 
be within the optimal region of the calorimeters with detector pseudorapidity |?7'^'^*| < 1.1 
or 1.5 < |?7'^'^*| < 2.5 |T^. When an electron is isolated, it is more likely to have originated 
from the decay of aW boson than from a b hadron. The efficiency of the combined electron 
identification requirements is ~ 60%. 



Jets, reconstructed with a cone algorithm of radius R = 0.5 [T^], must fail the electron 
requirements. The jet with highest transverse energy is required to have Et > 15 GeV and 
l^detj ^ 2 g_ rjj_^g second jet has to have Et > 10 GeV and \ri'^'^^\ < 4.0. Other jets in the 
event are counted if they have Et > 5 GeV and |?7'^'^*| < 4.0. We set the Et thresholds 
low and the \rj'^'^^\ region wide to maximize acceptance for signal; however, the efficiency to 
reconstruct jets close to the 5 GeV threshold is low. 

We identify a muon by the pattern of hits in the spectrometer drift tubes, and require 
an impact parameter < 20 cm between the spectrometer track and the primary vertex, a 
matching track in the calorimeter consistent with a minimum-ionizing particle, a matching 
central track, a signal in the scintillators surrounding the spectrometer within ±12 ns of the 
beam-crossing time, and penetration through one of the spectrometer toroids for momentum 
analysis. Most of these requirements are designed to reject cosmic rays and particles 
backscattered from the beamline magnets. Muons must be within the central region of 
the spectrometer, with \t]'^^^\ < 1.7. A muon is called "isolated" if Ai?(/i,jet) > 0.5 |15] for 
all jets with Et > 5 GeV. An isolated muon must have pt > 20 GeV and is attributed to 
the decay of a.W boson. A "tagging" muon has AR < 0.5 and > 4 GeV. It is attributed 
to the semileptonic decay of a 6 hadron in a jet, and thus identifies a b jet. The efficiency of 
the combined muon identification requirements is ~ 44% for isolated muons. 

Because a leptonically-decaying W boson is supposed to be present in each signal event, 
we require ^t > 15 GeV as evidence of a neutrino. 
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We use the NLO single top quark production cross sections to estimate that about 
66 s-channel and 153 t-channel events were produced at D0 during Run 1. Of these, we 
expect that about 15 s-channel and 35 t-channel events passed the trigger requirements and 
were recorded for analysis. 

Our analysis starts with a simple baseline selection of events that pass the triggers 
and have at least one reconstructed electron or isolated muon, and at least two jets with 
Et > 5 GcV and \r]'^'^*\ < 4.0. For the results presented in this Letter, we also demand at 
least one tagging muon ("//i"), to indicate the possible presence of a 6 jet. These minimal 
requirements reduce the ^ 1 million events in each channel to 116 e+jets//U events and 
110 /x+jcts//i events. The acceptance for single top quark events for these selections is 
0.2-0.3% per channel, which should yield 1 tagged event {tb and tqb, with electron and 
muon W decays combined). The expected number of events is small because the probability 
to identify at least one tagging muon in a single top quark event is only 6-11%. After these 
selections, 90% of the background in the electron channel is from QCD multijet production 
with a jet misidentified as an electron, 5% from ti events, and 5% from VF+jets (including 
1^1^ and WZ diboson events), where about two thirds of the VF+jets events have a light 
quark or gluon jet with a false tagging muon, and a quarter of the tagging muons are from 
c quark decays. In the muon channel, the background is 8% from W^+jets events, 6% from 
QCD bb events where a muon from a b decay mimics an isolated muon, and 4% from ti events. 
The remaining 82% of the background is from QCD multijet events with a coincident cosmic 
ray or beam-halo particle misidentified as an isolated muon. 

Next, we apply a set of loose criteria to remove mismeasured events and to reject 
backgrounds that do not have the same final-state characteristics as our signal. We reject 
events with more than one isolated lepton and any isolated photons. We remove events with 

close to 15 GcV and aligned with or opposite to a jet, or opposite an electron or isolated 
muon. We also reject events that have muons with clearly mismeasured pt- We require two, 
three, or four jets. To remove the remaining contamination from cosmic rays in the isolated 
muon channel, we reject events where the isolated muon and tagging muon are back-to-back; 
in particular, we require A0(isolyU, tag/x) < 2.4 rad. These criteria, together with the jet Et 
and I'q'^'^^l requirements and the threshold, reject 86% of the baseline multijet "electron" 
events, 95% of the cosmic ray and misreconstructed isolated "muon" events, 90% of the 
bb "isolated" muon events, 27% of the l^-|-jets events in the electron channel, 81% in the 
muon channel, and 55-73% of the ti background. The signal acceptances are reduced by 
14-51%. There remain 21 e+jets/jj, and 8 //-|-jets/// candidates in the data. 

Based on independent studies (see below), we apply the following requirements to obtain 
the best significance of signal over square-root of background in each channel: 

Electron Channel 

• ETiietl) + ET0et2) + Erie) + > 125 GeV 

• ^T(jet3) + 5 X Er(jet4) < 47 GeV 

• £;T(jetl) + 4 X > 155 GeV 
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Muon Channel 

• £^T(jetl) + ^r(jet2) + £;T(jet3) + EriietA) > 70 GeV 

• SrOetS) + 5 X ET(jet4) < 47 GeV 



The first criterion in each set was chosen by studying reconstructed compHEP [|1^] Monte 



Carlo (MC) W^+jets events, the second by examining HERWIG |]T3 ti MC events, and the 
third variable in the electron channel was determined from studies of QCD multijet data. 
The distributions were compared with signal MC from CompHEP. The cutoffs were optimized 
on combined samples of untagged and tagged events. Figure ^ shows the distribution of the 
second variable, designed to minimize tt background, for electron and muon events combined 
after all other selections have been applied. 
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FIG. 2. Variable used to reject tt background. 



After final selections, there is no evidence of an excess of signal over background, and we 
therefore use the results to set limits on the s-channel and t-channel single top quark cross 
sections. To do this, we must first determine the signal acceptance and the background in 
each channel. 

We obtain the signal acceptances using MC samples of s-channel and t-channel single 
top quark events from the CompHEP event generator, with the pythia package used 
to simulate fragmentation, initial-state and final-state radiation, the underlying event, 
and leptonic decays of the W boson. The MC events are processed through a detector 
simulation program based on the geant [|T^] package and a trigger simulation, and are then 
reconstructed. We apply all selections directly to the reconstructed MC events, except for 
several particle identification criteria, which we correct using factors measured in other D0 
data. Table | shows the acceptance for single top quark events after all selection requirements 
and corrections. 
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TABLE I. Signal acceptances (as percentages of the total 
cross sections), and numbers of events expected to remain 
after application of all selection criteria. 



Electron Channel Muon Channel 



Acceptances 

tb (0.255 ±0.022)% (0.112 ± 0.011)% 

tqb (0.168 ±0.015)% (0.083 ± 0.008)% 

Numbers of Events 

tb 0.18 ±0.03 0.08 ±0.01 

tqb 0.28 ±0.05 0.13 ±0.03 

W+iets 5.59 ±0.64 1.12 ±0.17 

QCD 5.92 ±0.58 0.40 ± 0.09 

tt 1.14 ±0.35 0.45 ±0.14 

Total Bkgd 12.65 ± 0.93 1.97 ±0.24 

Data 12 5 



The acceptance for ti pairs is calculated in a manner similar to that for signal and then 
converted to a number of events using the integrated luminosity for each channel and D0's 
value of the ti cross section 0. 

The QCD multijet background with a jet misidentified as an electron is measured using 
multijet data. The events are weighted by the probability that a jet mimics an electron 
for each jet that passes the electron Ej- and \ri'^'^^\ requirements. These probabilities are 
determined from the same multijet sample, but for < 15 GeV, and are found to be 
(0.0160 ±0.0016)% for < 1.1, and (0.0622 ± 0.0048)% for \r]'^''^\ > 1.5. We normalize 
the integrated luminosity of the multijet sample so as to match the data sample used in the 
search for the signal, and correct for a small difference in trigger efficiency between the two 
samples. 

The QCD bb background arises when both b quarks decay semileptonically to a muon, 
and one muon is misidentified as isolated. There are two ways for such events to mimic the 
signal. First, one of the b jets may not be reconstructed, and its muon can therefore appear 
to be isolated. Second, a muon can be emitted wide of its jet and be reconstructed as an 
isolated muon. The background from each source is measured using data collected with the 
same triggers as used for the muon signal. The events are required to pass all selections, 
except that the muon, which otherwise passes the isolated muon requirements, is within a 
jet. Events with truly isolated muons are excluded. Each event is then weighted by the 
probability that a nonisolated muon is reconstructed as an isolated one. This probability is 
measured using the same data sample, but for < 15 GeV, and found to be (2.94 ±0.53)% 
for the case of a "lost jet," and (1.38 ±0.25)% for a "wide /x," for muons with pt < 32 GeV. 
The probabilities are parametrized as a function of the muon pt- We calculate a weighted 
average of the two results to obtain the number of expected background events. 
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The background from W+jets is estimated by applying a set of tag-rate functions to 
untagged signal candidates that pass all final event selections. These tag-rate functions 
are measured using multijet data and correspond to the relative probability that a jet 
of given Et and 77*^^* has a tagging muon, for two run periods when the muon chambers 
had different operating efficiencies. We then correct the samples for a small difference in 
trigger efficiency between untagged and tagged events. We also correct the muon channel 
by a factor of 0.688 ± 0.034 to account for the effect of the A0 cutoff used to minimize 
cosmic ray backgrounds, a selection that cannot be applied directly. Finally, to avoid double 
counting, we subtract the fraction of events expected from tt and QCD backgrounds and 
single top quark signals. The remaining fraction of ly+jets in the untagged signal candidates 
is 66-92%, depending on the location of the electron or isolated muon. 

The numbers of events expected for the two signals and three backgrounds are shown 
in Table |, together with the final numbers of events in the candidate data samples, for the 
electron and muon channels. 

To calculate limits on the cross sections for single top quark production in the s-channel 
and t-channel modes, we use the numbers of observed events, the signal acceptances and 
backgrounds, and the integrated luminosities. Covariance matrices are used to describe 
the correlated uncertainties on these quantities. We use a Bayesian approach, with a flat 
prior for the single top quark cross section and a multivariate Gaussian prior for the other 
quantities. We calculate the likelihood functions in each decay channel and combine them 
to obtain the following 95% confidence level upper limits: 

• a{pp ^ tb + X) < 39 pb 

• a{pp tqh + X) < 58 pb 

To conclude, we have searched for electroweak production of single top quarks and find 
no evidence for such production. We set upper limits on the cross sections for s-channel 
production of th and t-channel production of tqh. The limits are consistent with expectations 
from the standard model. 
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